Physics of interface: Mott insulator barrier sandwiched between two metallic planes 
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We consider a heterostructure of a metal and a barrier with onsite correlation at half filling using 
unrestricted Hartree Fock. We find that above a certain value of correlation strength in the barrier 
planes, the system is a Mott insulator, while below this value the system still behaves like a gapless 
insulator. The energy spectrum is found to be very novel with the presence of multiple gaps. Thus 
the system remains non metallic for any finite value of correlation. 

PACS numbers: 73.20.-r, 71.27. +a, 71.30.+h, 71.10.Fd, 73.21.b, 73.40.C 
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INTRODUCTION 

Recently there has been a spurt in the study of multi- 
layered heterostructures, fabricated out of strongly cor- 
related materials. A wide range of systems have been 
studied experimentally and theoretically [3-0] • The inter- 
face of a band insulator and a Mott insulators can show 
metallic behavior [l| or can even become superconducting 
0. As shown by Thiel et al., such interfaces can be ma- 
nipulated by gate voltages thereby opening the prospect 
for interesting novel devices. 

In this letter, we investigate the interface of a metal 
and a system with finite onsite correlation. For onsite 
correlation greater than a certain value the system opens 
a gap, which mimics the Mott insulator. What hap- 
pens when we sandwich this system between two metallic 
planes? We find that the overall system behaves as an 
insulator not only in the region where there is a gap but 
also in the region where the system is gapless. 

The introduction of metallic planes induces multiple 
gaps in the energy spectrum for high values of onsite 
correlation U. This can be understood to be due to the 
emergence of vastly different types of sites both in plane 
and perpendicular to it. The gapless region behaves as a 
disordered insulator, with the disorder coming from the 
induced inhomogenity in the underlying Coulomb land- 
scape. There is a charge reconstruction across the planes 
for such a system as we vary U. In the region of param- 
eters where a gap opens up, this reconstruction is more 
vigorous. 

The interface that we have studied has been studied by 
others Is III using inhomogcneous dynamical mean field 
theory[8t-llll|. However, in the IDMFT approach used 
above all sites in each of the planes arc treated as idcnti- 
cal(spatial variation is taken only along the z direction), 
and thus misses out on crucial in plane modulations, 
which is one of the crucial reasons for the opening up 
of multiple gaps in the spectrum. It also considers a 
paramagnetic solution for the magnetic sector, which is 
not the correct solution for an underlying cubic lattice 



that we study. 

We employ the method of unrestricted Hartree Fock to 
solve the problem. It can capture the physics arising due 
to spatial variations both in and perpendicular to planes. 
It can handle the charge and magnetic sector on the same 
footing and treats correlations in a self consistent fash- 
ion. Recently we have shown the utility of this method 
in describing the metallic phase that arises in two dimen- 
sions as a result of competition between correlation and 
disorder (HEO!. 



Model and Method 

The barrier planes are described by the single orbital 
Hubbard model and the metallic planes by the non in- 
teracting tight-binding model. The Hamiltonian for the 
system is 
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Here the label a indexes the planes, and the label i in- 
dexes sites of the two-dimensional square lattice in each 
plane. The operator c\ a<7 (c iaa ) creates (destroys) an 
electron of spin a at site i on the plane a. We set the in- 
plane hopping t" to be nearest neighbor only, and equal 
to t, the hopping between planes, so that the lattice struc- 
ture is that of a simple cubic lattice. We take U a = U 
for the barrier planes, and zero for the metallic planes. 
The chemical potential // is calculated by demanding that 
there be exactly N electrons in the problem. This is done 
by taking the average of the N/2 th and the N/2 + 1 th 
energy level. We make no further assumption about the 
magnetic regime, and thus retain all spin indices in the 
formulas below. We label the barrier planes with a values 
from 1 to m. Thus a = l,m correspond to the metallic 
layers and all the other a values correspond to the barrier 
planes. 




FIG. 1: Top left and right: U = 4 and U = 5 for 28 x 28 x 5 
respectively. Middle left and right: U = 7 and U = 10 for 
28 x 28 x 5 respectively. Bottom left and right: U = 10 for 
28 x 28 x 4 and 28 x 28 x 6 respectively. 



FIG. 2: left: Gap at half- filling with increasing Mott Layers 
for (7 = 3 and 5 right :Effect of increasing number of layers on 
the multiple gaps. 



underlying nature of the overall system. 

For the quasi 2D geometry taken by us N = in x L 2 , 
where m is the total number of layers and L is the to- 
tal number of sites along both x and y direction, in the 
numerical work. We have performed finite size scaling of 
each of our results. We have shown the result of finite 
size scaling on the dc conductivity for a particular pa- 
rameter set(J7 = 5) in this paper due to paucity of space. 
The effect of varying the number of Mott layers has been 
studied in great detail in this paper. 



Calculated quantities 

The energy spectrum of the system is very rich, with 
the presence of multiple gaps in it. We have also cal- 
culated the charge profile, spin profile and the dc con- 
ductivity at zero temperature. We calculate the effect of 
correlations and layer variation on the multiple gaps. The 
value of charge and spin at a particular site in the cen- 
tral square of each plane is taken and is plotted against 
the plane index to clearly bring out the inhomogeneous 
profile along the z direction. The charge at a particular 
site is simply calculated as C = + rc^. The spin at 
a particular site is given by S — — ^1,4.) |- The dc 

conductivity is calculated using the Kubo formula, which 
at any temperature is given by: 

*M = 4 (2) 

N 6/3 ~ £a 

with A — ire 2 /hao, ao being the lattice spacing, and 
n a = Fermi function with energy e Q — /i. The f a p 
are matrix elements of the current operator j x = 
^Sit a( c l+xa o c i-<j ~ h.c), between exact single particle 
cigcnstates \ip a )> \ipp)> e ^ c > an d e a , ep are the correspond- 
ing eigenvalues. In this paper, conductivity/conductance 
is expressed in units of A = ire 2 /hao. 

We calculate the 'average' conductivity over small fre- 
quency intervals, Aw(Aw = rwj r , n = 1,2,3,4), and then 
differentiate the integrated conductivity to get cr{oS) at 
lj = nu r , n = 1.2.3[14j. We repeat the same calculation 
for each temperature slice upto kT — 0.1 The temper- 
ature dependence of the conductivity tells us about the 



Our Work and Results 

We have calculated the dc conductivity <7 zz (u; r ), within 
the accuracy allowed by these finite size systems. The 
conductivity is calculated in units of the universal con- 
ductance e 2 /h. Our method is able to capture the phe- 
nomenon of the gap opening as U goes above a certain 
threshold. We find two types of insulators, one is gapped 
and the other is not. This second type of insulating be- 
haviour which is shown for low U, U < U c , (U c is the 
value of U above which the system opens up a gap), is 
a disordered insulator, where the disorder is introduced 
by the introduction of the metal interfaces at the bottom 
and top. 

The effect of varying the insulating layer thickness on 
the charge order and spin order profile has been stud- 
ied in this work. As the gap at half filling closes the 
spin order parameter collapses to near zero values. The 
charge order parameter on the other hand varies less dra- 
matically across the gap closing transition. The charge 
profile builds up to its maximum value at the metallic 
planes and then dips to their lowest values in the first 
adjacent planes, from where there is massive charge de- 
pletion to the metallic planes, to gain Coulomb energy. 
As we move vertically deeper into the Mott planes, the 
total charge again picks up, reaching its bulk maxima at 
the central layer. There is overall symmetry around the 
central plane, as expected. As we increase the bulk thick- 
ness, we find that the charge peak at the central plane 
slowly approaches the value without the metallic planes, 
which is unity. This clearly indicates that the bulk be- 
comes gradually insensitive to the metallic planes at the 
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FIG. 3: Clockwise from top left:a Charge order variation vs 
layer index for a 18 x 18 x 10 system, for U = 2,3 and 4. b. 
Spin order parameter vs layer index for a 18 x 18 x 10 system 
for U = 2,3,4. c. Spin and charge order vs central site index 
for 17 = 4 as the system is squeezed along the z direction, d. 
Spin and charge order vs U for a 18 x 18 x 5 system. 



surfaces, as we increase the bulk layer thickness. 

Analysis of our results 

The sites having higher energy will have lower occupa- 
tion and vice versa. Fig la,b,c,d,e and f shows that differ- 
ent types of sites emerge in the problem that we study, as 
we systematically increase U. For different values of U, 
these gaps will open up successively. For m = 5, which 
means 3 intervening Mott layers, the first gap opens up 
roughly between U = 4 and U = 5, the second gap opens 
up close to U = 7 and the third gap opens up close to 
U = 10. The presence of 3 gaps clearly indicates the 
presence of broadly 4 different types of sites for U = 10. 

The 4 different types of sites correspond to the follow- 
ing. The first gap exactly at half filling corresponds to the 
situation when the U has become large enough to just in- 
troduce some anti ferromagnetic ordering(however weak) 
even in the metallic planes. Thus the overall system de- 
velops a rather inhomogeneous Neel ordering, where the 
inhomogeneity is along the z direction. Another way of 
looking at it is to say that in the absence of metallic 
planes, there would be a big gap at half filling. The in- 
troduction of metallic planes reduces this gap to a some- 
what lower value. As a result of this induction of anti 
ferromagnetic ordering(this will be discussed in more de- 
tails while explaining our spin/charge order results), two 
different types of sites emerge in the metallic plane for 
either spin. The reason for the generation of these two 
different type of sites in the metallic planes is as follows. 

To consider the emergence of two distinct types of sites 
due to the presence of interface, let us look at a particular 
spin(up/down). In the first type, an up electron sees 
comparatively high density of up spin and low density of 



down spin in the adjacent site in the Mott plane. Such an 
up electron will experience lower Hubbard repulsion and 
greater Pauli blocking if it wants to hop to the adjacent 
site in the Mott plane. An up electron in the second type 
of site sees comparatively lower up electron density and 
higher down spin density in the adjacent site on the Mott 
plane. This will lead to the second type of up electron 
experiencing higher Hubbard repulsion and lower Pauli 
blocking if it wants to gain kinetic energy by hopping to 
the adjacent Mott plane. 

As we increase U to about 7, the energy band of the 
electrons in the Mott planes decouple completely from 
the energy band of the electrons in the metallic planes, 
thus creating the second gap in the spectrum. The Mott 
layers now form the topmost band. As we crank up the 
U further, the environment seen by the sites in the Mott 
planes gets split further. An electron located in the cen- 
tral plane/s has the strongest Neel order, while an elec- 
tron located in the the Mott layers adjacent to the metal- 
lic planes have a relatively weaker Neel order. With in- 
creasing U the central layer becomes distinct from the 
two Mott layers which are adjacent to the metallic planes. 
This leads to the further splitting of the uppermost band 
at around U = 10. In Fig. le for system size 28x28x4, 
the uppermost gap is absent because of non existence of 
the central layer. 

Fig. 2a shows the plot of the gap at half filling as we 
increase the number of Mott layers for U = 3 and 5. It 
shows that for m = 4 and 5 there is no gap at half filling 
till U = 3. The gap opens up when U = 5. 

Fig. 2b shows how the multiple gaps in the energy spec- 
trum (A) behaves on increasing the number of Mott lay- 
ers for U = 10. This is because on increasing the width 
the number of sites at which the electron encounters U 
will rise, so the barrier becomes stiffer. 

Fig 3a shows the plots of total charge at a particular 
point in the central square of each layer vs the layer in- 
dex for m = 10, for three different values of U = 2,3,4. 
We can see clearly how as U increases the charge profile 
in the Mott layers become more and more non uniform 
along the z direction, thus developing a more pronounced 
hump. The charge depletion from the Mott layers just 
adjacent to the metallic planes and the the charge accu- 
mulation in the metallic planes increases with increasing 
U. In Fig 3b we have shown the plot of S vs layer index 
for m = 10 and U = 2,3,4. While S is almost zero and 
featureless for U — 2, it starts showing a prominent fea- 
ture which gets sharper as U is increased to 3 and then 
4. Fig 3c shows the plots of C and S for U = 4 but for 
three different values of m . — 5, 7, 9. As we increase the 
number of Mott layers, for a fixed U = 4, S becomes 
more pronounced, while C becomes more asymmetric in 
terms of increasing difference between the value in the 
central Mott plane and the Mott plane adjacent to the 
metal planes. 

In Fig 3d we show the plot of C and S vs U for 18 x 18 x 
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FIG. 4: Left: a vs kT for 30x30x5 system for U = 0,0.5,1,2. 
The qualitative difference in the curves as the system goes 
from metal to insulator with increasing U is captured very 
clearly. Right :Finite size scaling of <to vs kT for m = 5 and 
U = 5, for L=26,28 and 30. 



5 system. We have gone from U = to L 7 = 10 in integral 
steps. It clearly shows the emergence of three grossly 
different types of planes as we increase U. The three 
different types of planes are a) the two metallic planes, 
b) the two Mott planes adjacent to the metallic planes 
and c) the remaining Mott planc/s. We can see that as 
we increase U, more charge accumulates on the metallic 
planes and there is acute charge depletion from Mott 
planes adjacent to the metallic planes. There is charge 
depletion from central Mott plane/s also initially, but 
as we increase U further, the charge at the central Mott 
plane/planes increases slightly. This is accompanied by a 
large anti ferromagnetic ordering in the central planes/s. 
Thus the charges in the central Mott plane/s are unable 
to delocalize as it faces strong Coulomb repulsion on all 
sides and thus it piles up. 

In Fig. 4a we show the plots of dc conductivity against 
temperature T for U = 0,0.5,1 and 2 for a 30 x 30 x 5 
system, which is the largest system size that we have con- 
sidered. The u> r has been chosen to be twice the finite 
size spacing between energy levels for each of the system 
size that we have considered. The dc conductivity curve 
for C/ = 0, shows metallic behaviour as expected with a 
sharply falling profile with increasing T. We note that the 
value of conductivity obtained for the metal is severely 
suppressed. This is because, the highest contribution for 
the metal is at uj = 0, which cannot be sampled in our 
finite size simulation. Thus wc arc able to sample the 
<j(lo) at some low but finite u r where the conductivity 
has already fallen very sharply. As we increase U the dc 
conductivity for low T gets severely suppressed by sev- 
eral orders of magnitude and the system is an insulator 
even for U = 0.5. This is due to the emergence of sig- 
nificant amount of inhomogeneity in the charge density 
landscape at U = 0.5. As we increase U further, the 
amount of induced disorder increases further, suppress- 
ing conductivity even further. For U — 5, a small gap at 
half filling opens up, which seperates the occupied and 
unoccupied states. Fig 4b shows the finite size scaling 
effect on the c(w r ) vs T curves for U = 5, where we 
have taken L = 26,28,30, all of which show insulating 



behaviour. We have convinced ourselves that the result 
of finite size scaling on the gap at half filling, converges 
to a fixed value for L > 18 for m = 5. We find very 
high thermally activated conductivity which increases by 
several orders of magnitude. 



CONCLUSION 

A hetcrostructure system arising due to the sandwich 
of a barrier with finite width and onsite correlations be- 
tween two metallic planes is studied. There are multiple 
gaps in the spectrum, with increasing U . An inhomogc- 
neous spin and charge profile develops in the system. Dc 
conductivity calculations have been performed showing 
that the system is an insulator both in the gapped and 
gapless region. 
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